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We investigate the first-passage properties of bursty random walks on a finite one-dimensional 
interval of length L, in which unit-length steps to the left occur with probability close to one, while 
steps of length b to the right — "bursts" -— occur with small probability. This stochastic process 
provides a crude description of the early stages of virus spread in an organism after exposure. The 
interesting regime arises when b/L < 1, where the conditional exit time to reach L, corresponding 
to an infected state, has a non-monotonic dependence on initial position. Both the exit probability 
and the infection time exhibit complex dependences on the initial condition due to the interplay 
between the burst length and interval length. 



I. INTRODUCTION 



We are continually exposed to viruses. Despite these 
constant biological assaults, the immune system success- 
fully fends off most viruses. Considerable effort has been 
devoted to modeling the factors that influence whether a 
person exposed to a particular virus will eventually be- 
come ill jy . Typical theoretical models of viral infections 
account for the evolution of the number of infected cells, 
healthy cells, and viruses as a function of the rates of mi- 
croscopic infection and transmission rates. Such models 
have provided many useful insights about the dynamics 
of viral diseases H, 0] . 

In this work, we study a toy model — the bursty ran- 
dom walk (Fig. pi) — that captures one of the elements 
of viral infection dynamics. The position of the walk in 
one dimension represents the number of active viruses 
in an organism. Since the immune system constantly 
kills viruses, they are removed from the body at some 
specified rate, corresponding to steps to the left in the 
bursty random walk. However, with a small probability, 
a virus enters and successfully hijacks a cell, the outcome 
of which is a burst of a large number of new viruses into 
the host organism, corresponding to a long step to the 
right in the model. 

When the number of virus particles reaches zero, the 
organism may be viewed as being free of the disease. 
Conversely, when the number of viruses reaches a thresh- 
old value L, the organism can be viewed as either be- 
ing ill or dead. With this simplistic perspective, being 
cured or becoming ill is recast as a first-passage problem 
for the bursty random walk in an interval of length L. 
When the burst length b is small, the walk has a diffu- 
sive continuum limit whose first-passage properties are 
well known p, [1 . However, if the burst length is of the 
order of the system length, this burstiness effects strongly 
affect the first-passage characteristics. This large-burst 
limit should be applicable to infectious processes where 



the threshold number of viruses for being ill is not large 
and the number of new viruses created in a burst event is 
a finite fraction of this threshold || . Related discreteness 
effects were found in the first-passage characteristics of a 
random walk that hops uniformly within a range [—a, a] 
in the interval [0, L], with a < L M. 




FIG. 1: A bursty random walk with burst length 6 = 6. 

In the next section, we define the model and the ba- 
sic fi rst- passage quantities that we will investigate. In 
Sees. Ill and [V, we determine the exit probabilities and 



the average exit times to either end of the interval as a 
function of the burst length b. When b < L, very different 
first-passage properties arise compared to those for pure 
diffusion in the interval. Perhaps the most striking is the 
conditional exit time to reach x = L, corresponding to 
a state of infection, which has a non-monotonic depen- 
dence on the starting position x. We compute these first- 
passage properties from the backward Kolmogorov equa- 
tions for the exit probabilities and exit times p, m . In the 
concluding section, we briefly discuss the corresponding 
first-passage properties for the bursty birth/death model. 
This process accounts for the feature that bursts should 
occur at a rate that is proportional to the number of live 
viruses. It is naturally to model this situation by defining 
the rate at which steps occur to be proportional to the 
current position of a bursty walk on the interval. 



II. THE MODEL 



In the bursty random walk, unit-length steps to the 
left occur with probability q, while long steps (bursts) of 



length b occur with probability p = 1 — q (Fig. [l]). We 
choose p and q so that the average position of the walk 
does not change at each step; however, most of our re- 
sults are derived for general p and q. The motivation for 
considering these hopping probabilities is based on the 
experimental observation that viral counts in an organ- 
ism often remain nearly constant for time periods much 
longer than the lifetime of individual viruses. Such a near 
constancy could only arise if an organism produces new 
viruses (by bursts) and clears viruses at similar overall 
rates gj. 

With the constraint that the number of virus particles 
remains fixed, on average, the respective probabilities of 
making a single step to the right and to the left are 

1 b 
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The bursty random walk is confined to the finite inter- 
val [0,L], where the coordinate represents the number 
of live viruses. The state where the virus is cleared is 
represented by the point x = 0, while the state where a 
sufficient number of viruses exists that the organism is ill 
or dead is represented by x — L. Our goal is to under- 
stand first-passage properties of this bursty walk that are 
relevant to the state of health of the organism. Namely, 
what is the probability that the organism becomes cured 
or becomes ill, corresponding to the walk eventually exit- 
ing through the left edge or the right edge of the interval, 
respectively? What is the time needed for the organism 
to become cured or ill? 

To set the stage for our results, let us recall some well- 
known first-passage properties for the isotropic nearest- 
neighbor random walk on the interval 0. Let £+(x) 
denote the probability that this random walk eventually 
exits the interval at L without ever touching x = 0, given 
that the walk started at an arbitrary point x, with < 
x < L. The complementary exit probability to the left 
boundary is £-{x) = 1 — £+(x). These exit probabilities 
satisfy the recursion 



£±{x) = \£±{x - 1) + \£±{x + 1) 



(2) 



subject to the boundary conditions £+(0) = 0, £+(L) = 
1, or £_(0) = 1, £-{L) = 0. This recursion expresses 
the exit probability starting at x as the probability of 
first taking a step to the left or right (the factor 1/2) and 
then exiting from either x — 1 or x + 1, respectively. The 
solution to Eq. (0) with these boundary conditions is: 



subject to the boundary conditions i(0) = t(L) = 0. 
Equation (||) expresses the average exit time from x as the 
time for the first step (the additive factor 1) plus the exit 
time from the new positions (either x =b 1); the factor i 
accounts for the probability for each of these two choices. 
Similarly, we also define the conditional exit times, t±(x), 
as the average times for the walk to leave the interval by 
the right or the left boundary, respectively, without ever 
reaching the opposite boundary. These conditional exit 
times satisfy [|] 

C±{x) = \C±(x-l) + \C±(x+l) + £±{x), (5) 

with C± = £±t±, and this equation is subject to the 
boundary conditions C±(0) = C±(L) = 0. For the 
nearest-neighbor random walk, the exit times are given 
by 



t(x) = \ x(L — x) , 
t + (x) = l(L 2 -x 2 ), 
t_( x ) = \{2Lx-x 2 ). 



(6) 
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Our goal is to determine the results analogous to 
Eqs. (pi) and (0) for the bursty random walk. As we 
shall see, first-passage properties depend only on bjL as 
long as this ratio is nonzero. 



III. EXIT PROBABILITIES 



For the bursty random walk with burst length b, we 
may naturally define two distinct types of exit probabil- 
ities to the right boundary: 



• the total exit probability £+ (x) that the walk even- 
tually reaches any point at the right boundary or 
beyond, without ever touching the left boundary 
x = 0] 

• the restricted exit probability lZ m (x) that the walk 
eventually reaches the specific point L + m (with 
< m < b — 1), without ever touching the left 
boundary or any other point beyond the right 
boundary. There are b such distinct restricted exit 
probabilities, lZ m (x), with m = 0, 1, . . . , b — 1. 



We also define t(x) as the average time for the walk to 
leave the interval at either end when it starts at x. This 
unconditional exit time satisfies the recursion 



t(x) = | t(x - 1) + I t(x + 1) + 1 . 



(4) 



While the total exit probability is most relevant physi- 
cally, because it corresponds to the probability of illness 
for a given level of initial exposure, the restricted exit 
probabilities display intriguing features that stem from 
the bursty character of the walk 1t| . 



The total exit probability to the right boundary sat- 
isfies the recursion 



tion aX x + bx + c, we obtain 



£+{x) = q£+{x - 1) + p£ + (x + 6) , 



(7) 



that represents the extension of Eq. (J2j) to the bursty 
random walk. This recursion expresses exit via the right 
boundary, when starting from x, either by taking the 
first step to the left (probability q) , after which exit from 
x — 1 occurs, or by first stepping to the right (probability 
p), after which exit from x + b may occur. This recur- 
sion must be supplemented by the boundary conditions 
£+(x) = 1 for all x > L and £+{0) = 0. Namely, a walk 
that starts at x > L has already exited, while a walk 
that starts at x = can never exit via the right bound- 
ary. The equations for the restricted exit probabilities 
are similar to d7h, but are now subject to the boundary 
conditions 7?. m (0) = and lZ m (L + k) = Sk, m - 

While the exit probabilities can be obtained by enu- 
merating all random walk trajectories to the exit point 
and computing the probabilities for all these paths, the 
above recursions provide the same results much more eas- 
ily [p, El. We will use different methods to solve Eqs. (M) 
for short and large burst lengths, and therefore study 
these cases separately. 



A. Burst lengths 6 = 2,3, 



For the first non-trivial case of burst length b — 2, we 
solve the constant-coefficient recursion (Q) by attempting 
a solution of the form E + (x) = X x . This leads to the 
characteristic equation A 3 — 3A+2 = 0, with solutions A = 
—2 and A = 1 (doubly degenerate). Henceforth, we use A 
to denote the first root of the characteristic polynomial. 
The general solution to Eq. (0) thus is £+{x) = aX x + 
bx + c. Invoking the boundary conditions, the total exit 
probability to the right boundary is 
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(A l h 



zpw 



for the restricted exit probabilities to x — L and to 
x = L + l, respectively. Parenthetically, once we 
know one of IZo(x) or IZi(x), the other is determined 
by the martingale property that the mean position of 
the bursty walk always remains fixed [|3|. That is, af- 
ter all the probability has reached an absorbing bound- 
ary, the two restricted exit probabilities are related by 
Ox[l-TZ (x)-'R,-i(x)}+Lxno(x) + (L+l)xn 1 (x) = x. 
The restricted exit probabilities initially grow nearly lin- 
early in x/L (Fig. 0), but then oscillate violently as 
x — > L. The total exit probability £+{x) is a nearly 
linear function of x for small x but its slope develops 
oscillations as x — > L. 

This same calculational method can be extended to 
longer bursts. By assuming an exponential solution of the 
form £+{x) — \ x in Eq. (0), the characteristic polynomial 
generically is (A— 1) 2 A(X), where A(X) is a polynomial of 
order b — 1. Explicit closed-form solutions can therefore 
be obtained for b < 5, but numerically exact results can 
be obtained for any burst length; details for the case 
b = 3 are given in appendix KL Typical results are shown 
in Fig. H for burst lengths b — 2, 10, and also 6 = 4, 
4, and L. The total exit probability is very close to 
linear function with slope less than one for x < L — b, 
but deviates from linearity within one burst length from 
x = L. The restricted exit probabilities are also nearly 
linear functions for x < L — b, but oscillate violently in 
the boundary region. 
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(8) 



B. Long Bursts 



For the restricted exit probabilities to a specific point 
in the range [L, L + b — 1], the boundary conditions are: 

K o (0)=K o (L + l) = Q, Tl (L) = l, exit to x = L, 
TZ 1 (0)=TZ 1 (L) = 0, TLi{L + l) = l exit to x = L+l. 

Applying these boundary conditions to the general solu- 



When the burst length is of the order of the inter- 
val length, we can simplify the determination of the exit 
probabilities by considering separate recursions in each 
of the disjoint subintervals [L — b,L), [L — 26, L — b — 1], 
[L — 36, L — 26— 1], etc., instead of directly solving for the 
roots of a characteristic polynomial of order 6 — 1 . As we 
shall see, this partitioning significantly reduces the order 
of the recursions for the exit probabilities. 



1. Total exit probabilities 



In the extreme situation where the burst length b > L, 
a single burst results in exit at or beyond the right end of 
the interval. Thus the total exit probability satisfies the 
recursion £+(x) — q£ + (x — 1) + p; that is, either the walk 
steps to the left and then exits from x — 1, or the walk 
steps to the right and exits immediately. The solution to 
this recursion is a constant plus an exponential function. 
The boundary condition £+(0) = immediately gives 



£ + (x) = l-q a 



(10) 



(11) 



Because of the overwhelming probability of stepping to 
the left, the exit probability to the right boundary is not 
close to one for x — >• L from below. As an example, for 
b = L, wehave£+(L-l)^ 1-e" 1 w 0.6321 (Fig.|(c)). 

For the case L/2 < b < L, we partition [0,L] into 
the subintervals [0, L — b — 1] (defined as region I) and 
[L — b,L] (region II), as indicated in Fig. |[ Making a 
slight abuse of notation, we define £ l (x) and £ Il (x) as 
the total exit probabilities to x > L, when starting at a 
point x that is in either region I or region II, respectively. 
These exit probabilities satisfy 

£ l {x) = q£ I (x-l)+p£ II (x + b), 

£ u {x) =q£ u (x-l)+ P . 

These recursions are identical in form to Eqs. (0), but 
with the subinterval explicitly identified. Thus, from 
example, exit to x > L, when starting from a point 
x in region I, can occur by taking a step to the left 
with probability q and then exiting from x — 1 (neces- 
sarily in region I), or by taking a step to the right with 
probability p and then exiting from x + b (necessarily 
in region II). Equations (O) are subject to the bound- 
ary condition £ l (0) = as well as the joining condition 
£ il {L-b)=q£ l {L-b-l)+p. 

By this partitioning, the exit probabilities in each 
subinterval are functionally distinct and can be solved 
separately. In the second of Eqs. (pd|), a particular solu- 
tion is £p al — 1. Thus the general solution has the form 
£ u (x) = 1 + Aq x . Substituting this expression in the 
first of Eqs. ([ll]), now gives the closed recursion £ l {x) — 
q£ l (x — 1) + p + Apq x . With the inhomogeneous term 
p + Apq x , the general solution is £ l (x) = A+(Bx + C)q x . 
Using the boundary condition £ l (0) — 0, and substitut- 
ing this form for £ l {x) into the first of Eqs. dO), we find 
B = Apq b . Finally, we invoke the joining condition and 
obtain 



£\x) = l-q 1 



xpq 



x+b 



£ u (x) = 1- 



, (12) 



i-ypq b ' l-ypq b 

where y = L—b—1. Notice again that because of the large 
probability of hopping to the left, £\\{x) is discontinuous 
as x — > L. 
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FIG. 2: Exit probabilities for the bursty random walk for: (a) 
burst length b = 3, (b) 6 = 10, and (c) b = L, L/2, and L/2,. 
Simulations are on a system of length L — 100. 



For L/3 < b < L/2, we partition [0,L] into the three 
subintervals [0,1-26-1], [L-2b,L-b-l],aod [L-b,L], 
(regions I, II, and III respectively) and solve the general- 
ization of Eqs. (11) to three intervals, supplemented by 
two joiningconditions at x = L — b and at x = L — 26 
(appendix |BJ) . As shown in Fig. g, the total exit proba- 
bility has two (barely visible) singularities and deviates 
considerably from linearity within one burst length from 
the right boundary. Generally, for a partitioning into k 
intervals, the slope of the total exit probability is discon- 
tinuous at the boundary between intervals k and k — 1, 
the second derivative is discontinuous at the boundary 
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FIG. 3: Partitioning of the interval into [0, L — b — 1] (region 
I) and [L — b, L] (region II). To exit from region I requires at 
least two bursts. 



between intervals k — 1 and k — 2, the third derivative 
is discontinuous at the boundary between intervals k — 2 
and k — 3, etc. A similar intricate pattern of a sequence of 
progressively weaker singularities arises in various frag- 
mentation models [pi [HI . 
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FIG. 4: Restricted exit probabilities to x — L for the cases of 
b = L/2 and b — L/3. The arbitrary vertical scale has been 
fixed by setting the integral of 7?.(a;) over the interval equal 
to 1. 



2. Restricted exit probabilities 



The restricted exit probabilities to a specific point un- 
dergo a more dramatic sequence of discontinuities be- 
tween successive subintervals. We again start with the 
case where b lies in the range [L/2,L] so that there are 
two subintervals to consider: [0, L — b — 1] and [L — b, L]. 
For concreteness we determine the exit probability to the 
specific site x = L; similar behavior arises for other exit 
points in [L, L + b — 1]. Now the recursion relations for 
the restricted exit probabilities are 



TZ\x) = qTZ'ix - 1) + pTZ u (x 
K l \x) =qU u {x-l). 



b). 



(13) 



Since we seek only the exit probability to i = L, we 
simplify notation by omitting the subscript that specifies 
the exit location; thus TvLo — > 1Z. The first equation states 
that to reach x = L from subinterval I, the walk can 
either step left (probability q) and exits from x — 1, or 
the walk steps to the right (probability p) and exits from 
x + b. The second equation states that to reach x = L 
from within subinterval II, the only possibility is to step 
to the left; a burst would lead to exit at a point x > L, 
which does do not contribute to the exit probability to 
x = L. The recursions (O) must be supplemented by 
the boundary condition 1Zr = and the joining condition 
K U {L -b) = qTZ l (L-b-l)+p. Notice that exit to L 
can occur only if the walk is at the point x = L — b. 



at the boundary between intervals k and k — 1, the first 
derivative is discontinuous at the boundary between in- 
tervals k — 1 and k — 2, etc.; the pattern is similar to that 
for the total exit probability, but the discontinuities are 
more prominent here since they begin with the function 
itself rather than with the first derivative. 



IV. FIRST-PASSAGE TIMES 



By adapting Eq. (Eh to the bursty random walk, the 
unconditional mean first-passage time satisfies 



t{x) = qt(x-l) +pt(x + b) + 1 



(15) 



subject to the boundary conditions t(0) — and also 
t(L + m) = for m = 0, 1, . . .b — 1. Similarly, the quan- 
tities C±(x), which are related to the conditional exit 
times, satisfy the recursion (see Eq. (jj|)) 

Similarly, the conditional mean first-passage times 
C±(x) satisfy the recursion (see Eq. (||)) 



C±(x) = qC±(x— 1) + pC±(x+b) + £±(x) , 



(16) 



subject to the same boundary conditions as for t(x) itself. 
Again, we treat the exit times separately for short and 
for long bursts. 



Employing the same method as that used to obtain 
Eqs. (|l2|), we now obtain obtain 
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(14) 



This solution method can be extended to more subinter- 
vals; the results for b = L/2 (two intervals) and b = L/3 
(three intervals) are shown in Fig. 0. For a partition- 
ing into k intervals, the exit probability is discontinuous 



Short bursts 



For the first non-trivial case of b — 2, let us focus on 
the unbiased case of p = | and q = | for simplicity. We 
solve the recursion (UB) with these values of p and q by 
noting that the inhomogeneous term can be eliminated 
by writing t(x) = T(x) — x 2 /2. Substituting this ansatz 



into Eq. (jig), we find that T{x) obeys this same equation, 
but without the inhomogeneous term. From our analysis 
of the exit probability in Sec. Ill A[ the general solution 
is T{x) = a\ x + bx + c, with A = —2, subject to the 
boundary conditions T(0) = 0, T(L) = L 2 /2, T(L + 1) = 
(L + l) 2 /2 that correspond to i(0) = t(L) = t(L + 1) = 
0. We thereby obtain, for the unconditional mean first- 
passage time, 



*(«) = 2 X ( L - X )- 



l {L + l)[{\*-l)-l{\ L -l) 
2 (A i + 1 -l)-^±i(A L -l) 



(17) 



with A = — 2. The second term represents a tiny correc- 
tion to the leading diffusive behavior of \x(L — x). 



B. Long Bursts 



In the extreme case of burst length b > L, the walk 
exits after any single burst, and the unconditional first- 
passage time satisfies t(x) = qt(x — 1) + 1, subject to the 
boundary condition t(0) = 0. The solution is 



t{x) = 



l-q x 
1-5 



(18) 



Similarly, the conditional mean first-passage time to the 
right boundary, t + = C + /£ + , is determined from the re- 
cursion 

C+{x) = qC+{x-l)+£ + {x) = qC+{x-l) + l-q x , (19) 

subject to the boundary condition C + (0) = 0. The solu- 
tion now is 



t+(x) 



1 



ig 



1 — q 1 — q x 



(20) 



The conditional exit time £_ may be obtained from the 
conservation statement t(x) = £~(x)t-(x) + £ + (x)t + (x) 
and gives t-(x) = x. An apparently paradoxical fea- 
ture is that the exit time t + increases when the starting 
point is closer to x — L (Fig. ||(d)). This behavior arises 
because steps to the left occur with overwhelming prob- 
ability. Thus a walk that starts near x — L will almost 
surely hop a considerable distance to the left before a 
burst occurs. However a walk that starts near x = can 
only hop a short distance to the left before a burst must 
occur to ensure exit at the right boundary. 

For L/2 < b < L, we again partition the interval into 
the subintervals [0, L — b — 1] (region I) and [L — b, L] (re- 
gion II) and denote the mean first-passage times within 
each as t l (x) and t (x) respectively. The unconditional 
mean first-passage time satisfies 



t l {x) = qt 1 (x - 1) + pt u (x + b) + 1 
t u ( x ) = qt u (x - 1) + 1 , 



(21) 
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FIG. 5: Normalized mean first-passage times r = t/(L 2 /D), 
with D = b/2, for: (a) the nearest-neighbor random walk, 
and the bursty random walk with: (b) burst length b — L/5, 
(c) b = L/2, and (d) b = 2L. Shown are the unconditional 
first-passage time t(x) and the conditional times, t±(x), to 
the right and left boundary, respectively. 



subject to the boundary condition t l (x) = and the join- 
ing condition t ll (L — b) = qt l (L — b — 1) + 1. Solving first 
for i 11 and then using this solution in the equation for t 1 , 
we obtain 



t\x) = 
t u (x) = 



x(q~ 



,6+x 



2(f 



1 - ypq b 
[1 - ypq b + 



2q x ) 



v (i - ypq b ) 



with y = L — b — 1 . 



(22) 



decreases as x increases. The source of this decreasing 
dependence on x in this range is that a particle with 
x > L — b is increasingly likely to reach a point that is 
less than L — b as x decreases toward L — b. Once the 
point x = L — b is crossed, two bursts are required for 
exit to the right and typically there will be many steps 
to the left between these two bursts. Thus the exit time 
increases rapidly as the starting point approaches L — b 
from above. 



For the conditional first-passage time to the right 
boundary, the quantity C+(x) = £+(x)£+(x) satisfies 



C\x) = qC\x - 1) + pC u (x + b)+ £\x) , 
C u (x)=qC I (x-l) + £ u {x), 



(23) 



subject to the boundary condition C l (x) = and the 
joining condition C U (L - b) = qC\L - b - 1) + S U (L - b). 
Again, we have made the notational abuse of dropping 
the subscript ± and focusing only on the exit time to the 
right boundary. Solving these equations for C+ and di- 
viding by £_|_ (x) yields the conditional first-passage time 
to the right boundary (Fig. |1). This same calculation can 
be straightforwardly (but tediously) extended to smaller 
values of b, corresponding to more subintervals. 




FIG. 6: The normalized conditional exit time to the right 
boundary, r+ = t+/(L 2 /D) for the case b/L = 0.8. 

A peculiar feature of the conditional first-passage time 
t+{x) is its non-monotonic dependence on x as the burst 
length becomes of the order of the system length (Fig. H). 
This non-monotonicity has a simple origin. For x/L < 
1, a typical walk will move a considerable distance to 
the left before exit occurs. Thus, in some sense, points 
near the right boundary are "further" from the exit than 
points in the interior of the interval. Similarly, a particle 
that starts near x — must quickly hop to the right 
to avoid exiting at the left boundary. Thus again, the 
exit time to the right is an increasing function of x in 
this range. Finally, for a particle that starts in a narrow 
range in which x is slightly larger than L—b, the exit time 



DISCUSSION 



We investigated the first-passage properties of the 
bursty random walk on a finite interval, where short steps 
to the left occur with a high probability, while long steps 
to the right — "bursts" — occur with a small probabil- 
ity. The disparity in these hopping probabilities is needed 
to ensure that there is no net displacement of a random 
walker, a feature that maximizes the time for the walker 
to survive within the interval. This model was motivated 
by the problem of the early stages of virus spread after 
initial exposure B. 

When the burst length is short, there are only small 
corrections to the well-known first-passage properties of 
the nearest neighbor random walk. Conversely, when the 
burst length is of the order of the interval length, dis- 
creteness effects play an important role. For such burst 
lengths, we solved for first-passage properties by parti- 
tioning the full interval into disjoint subintervals of length 
b, solving each one separately, and then patching together 
these subinterval solutions by invoking appropriate join- 
ing conditions. Strikingly, the mean first-passage time to 
the right boundary, corresponding to the time for a host 
organism to become ill, has a non-monotonic dependence 
on the initial location for b/L < 1 (Fig. H). Another basic 
feature of the first-passage properties for large b is that 
they are functions of b/L rather then depending on b and 
L separately. 

In spite of the strange behavior of the mean first- 
passage time, the distribution of first-passage times is 
generically characterized by an exponential decay, but 
with superimposed oscillations due to burstiness. Con- 
sequently, higher moments of the first-passage times can 
be simply characterized by powers of the first moment. 

If one takes seriously the equivalence that the position 
of the random walker as equivalent to the number of ac- 
tive viruses, then the frequency of bursts as well as the 
frequency of virus death events should also be propor- 
tional to the current position of the walk. Thus it would 
be more realistic to consider the bursty birth/death pro- 
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FIG. 7: The normalized unconditional and conditional exit 
times, r = t/L, for the bursty birth/death process with b/L = 
0.75. 



cess, where the rate at which the random walker hops is 
proportional to its current location. If a step does occur, 
then a unit-length step to the left occurs with proba- 
bility q and a step of length b to the right occurs with 
probability p <C q. 

Because the exit probabilities are independent of the 



rate at which steps occur, all our results about exit prob- 
abilities continue to hold for the bursty birth/death pro- 
cess. However, exit times for bursty birth/death are quite 
different from those of the bursty random walk. For ex- 
ample, the unconditional exit time for bursty birth/death 
satisfies the recursion 



t(x) = qt{x - 1) + qt(x + b) + St(x) , 



(24) 



where St(x), the microscopic time step at position x, is 
proportional to 1/x. For the classic birth/death pro- 
cess (burst length b = 1) and in the continuum limit, 
Eq. ( pil ) becomes t"(x) = — 2/x with solution t(x) = 
2xla(L/x). Over most of the interval range, this exit 
time scales linearly with L, compared to t(x) ~ L 2 for 
the exit time of the nearest-neighbor random walk. For 
bursty birth/death, representative results for exit times 
are given in Fig. [7|. While no longer non-monotonic in 
x, the conditional exit time t+(x) has a near plateau 
when the initial position x < L — b and then decreases 
in x. Thus once an infection has progressed to a certain 
threshold, illness quickly ensues. 
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Appendix A: Exit probabilities for burst length 6 = 3 



For burst length b = 3, the recursion relation for the total exit probability to the right boundary is 



(Al) 



Assuming the exponential form £ + = X x and substituting into (Al), the characteristic equation is A 4 — 4A + 3 = 0, 
with solutions A = 1 (doubly degenerate) and A± = — l±i^/2 = y/Se 1 ^, with <fi = tan -1 (—a/2)- The general solution is 
£+{x) = aX++b\ x _+cx+d. Now we impose the boundary conditions £ + (0) =0and£ + (L) = £+(L+l) = £ + {L+2) = 1 
one by one. The boundary condition £+(0) = gives 



£+{x) = a{\ x + - 1) + b{\ x _ - 1) + ex . 
The boundary condition £+{L) — 1 gives 



£ + (x) 



(XI 1) - (X L + 1) 



aa(x) -f ba*(x) + — , 



L 



+ b 



(x x _ - 1) - (A 1 : - 1)- 



x 
L 



with a(L) = 0. Next we impose £+(L + 1) = 1 to give 



£ + {x) = a[a(x)a*(L + 1) - a*(x)a(L + 1)] + g + (l - £) ^+1) 

x ( x\ a*(x) 
^ aW (x,L + l) + I + (l- I )^ Yy 



(A2) 



(A3) 

(A4) 
(A5) 



where the Wronskian W(L + 1,L + 1) = 0. Finally, imposing £+(L + 2) = 1 gives 



£+(x) 



W(x,L + l) 
W{L + 2,L + 1) 



L + 2 la*(L + 2) 



La*(L + l) 



1 a*(x) 



L La*(L + l) 



By inspection, it is clear that Eq. (|A6[) satisfies all the boundary conditions; this solution is also real. 



(A6) 



A similar calculation can be performed for all the restricted exit probabilities. For example, for the restricted exit 
probability 1Zq{x) to L, we start with Eq. (|A2j) and next impose TZq(L + 1) = to give 



lZo(x) = a 



(A* - 1) - (A^ +1 - 1) 



L + l 



(XL - 1) - (XL +i - 1) 



aP(x) + b/3*{x), 



(A7) 



with /3(L + 1) = 0. The boundary condition TZq(L + 2) = leads to IZo(x) = aV(x,L + 2), where the Wronskian is 
now defined as V(x,y) = f3(x)f3*(y) — /3*(x)/3(y). Imposing the boundary condition TZq(L) = 1 gives the final result 



n {x) 



V(x,L + 2) 



V(L,L + 2) ' 
For the other two restricted exit probabilities, the same calculation as that outlined above gives 

W{x,L + 2) ^ ,_^ W(x,L + l) 



K x {x) 



n 2 {x) 



W{L+l,L + 2) ' zv ; W(L + 2,L + l) 

These results for the total and restricted exit probabilities are plotted in Fig. 0. 



(A8) 



(A9) 



Appendix B: Exit probabilities for burst length L/3 < b < L/2 



When the burst length b is in the range [^, -j], the interval naturally divides into the three subintervals [0, L — 26— 1], 
[L — 2b, L — b — 1], and [L — b, L]. The recursion relations satisfied by the total exit probability to the right edge of 
the interval are: 

£\x) = q£\x - 1) + P £ u (x + b) , 

£ l \x) = q£ u {x - 1) +p£ m (x + b) , (Bl) 

£ ll \x)=q£ ll \x-l)+p. 

These exit probabilities must also satisfy the joining and boundary conditions 

£ l (0) = 0, 
£ u (L-2b) = q£ I {L-2b-l)+p£ m (L-b), 
£ m {L-b) = q£ u (L-b-l)+p. 



We generalize the approach used to solve the two-interval case (cf. Eq. (|T|)) by first solving for £ IU in the form 
£ lu = 1 + Aq x , substituting this result into the recursion for £ u to obtain its general form, and finally substituting 
the result for £ u into the recursion for £ l . All the unknown constants may then be fixed by the boundary and joining 
conditions, and the final result is: 

= x _ q x {pq b [2(x -y)+ pq b (b + x - y)(b + X - y + 1)] + 2} 
pq b [(b~y)(b-y+l)pq b -2y]+2 
£ U {x) = l 2q*[pq b (x-y) + l] 



£ lll (x) = 1 



pq b [(b -y)(b-y+ l)pq b - 2y] + 2 ' 

2q x 
pq b [(b-y)(b-y + l)pq b -2y] +2 ' 



(B2) 
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where y = L — b — 1. This procedure can be continued to as many subintervals as desired both for the total and for 

the restricted exit probabilities. 
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